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We suggest a constrained instanton (CI) solution in the physical QCD vacuum which is described 
by large-scale vacuum field fluctuations. This solution decays exponentially at large distances. It 
is stable only if the interaction of the instanton with the background vacuum field is small and 
additional constraints are introduced. The CI solution is explicitly constructed in the ansatz form, 
and the two-point vacuum correlator of gluon field strengths is calculated in the framework of the 
effective instanton vacuum model. At small distances the results are qualitatively similar to the 
single instanton case, in particular, the form factor Di is small, which is in agreement with the 
lattice calculations. 
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I. INTRODUCTION 



The non-perturbative vacuum of QCD is densely populated by long - wave fluctuations of gluon and quark fields. 
The order parameters of this complicated state are characterized by the vacuum matrix elements of various singlet 
combinations of quark and gluon fields, condensates: (: qq :), (: F"F® :),(: q{a tlv F^ v ^-)q :) , etc. The nonzero 
quark condensate (: qq :) is responsible for the spontaneous breakdown of chiral symmetry, and its value was estimated 
a long time ago within the current algebra approach. The nonzero gluon condensate (: F^F® :) through trace 
anomaly provides the mass scale for hadrons, and its value was estimated within the QCD sum rule (SR) approach. 
The importance of the QCD vacuum properties for hadron phenomenology has been established by Shifman, Vainshtcin 
and Zakharov jlj. They used the operator product expansion to relate the behavior of hadron current correlation 
functions at short distances to a small set of condensates. The values of low - dimensional condensates were obtained 
phcnomcnologically from the QCD SR analysis of various hadron channels. 

Later on the nonlocal vacuum condensates or vacuum correlators have been introduced They describe the 

distribution of quarks and gluons in the non-perturbative vacuum. Physically, it means that vacuum quarks and 
gluons can flow through the vacuum with nonzero momentum. From this point of view the standard vacuum 
expectation values (VEVs) like (: qq :), (: qD 2 q :) , (: g 2 F 2 :) , ... appear as expansion coefficients of the quark 

M(x) = ^: q(0)E(0, x)q(x) :\ and gluon D pv - p<J {x) correlators in a Taylor series in the variable x 2 /4. 

The correlator D fiu ' pa ' (x) of gluonic strengths^], 

£>~(x - y) = (: TrF^(x)E(x, y)F p °(y)E(y, x) :) , (1) 

may be parameterized in the form consistent with general requirements of the gauge and Lorentz symmetries as Q : 

D^(x) = ^(-F 2 :) {(<W<W - 5^5 vp )[D(x 2 ) + D^x 2 )] + (2) 

1 K^p^pOva 'u^v p 1 X V X (jO Xi/XpO^fjJ Qx 2 



1 We follow the convention when the coupling constant is absorbed into the gauge field A v (x). 
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where E(x,y) = Pexp(i J A ll (z)dz fl ) is the path-ordered Schwinger phase factor (the integration is performed 

A a A a 

along the straight line) required for gauge invariance and A^(z) — A^(z) — , F^(x) = F° v (x) — , F£ v (x) = 

d/j,A*(x) — d v A"(x) + f abc A h 1 (x)A c v (x). The P— exponential ensures the parallel transport of color from one point 
to other. In (^), (: F 2 :) = (: F°(0)F°(0) :) is a gluon condensate, and D(x 2 ) and Di(x 2 ) are form factors which 
characterize nonlocal properties of the condensate in different directions. They are normalized at zero by the conditions 
D(0) = k, -Di(O) = 1 — k, that depend on the dynamics considered. For example, for the self-dual fields k = 1, while 
in the Abelian theory without monopoles the Bianchi identity provides k — 0. 

The gluon correlators D^ l/pa (x) are involved in an analysis M of the spectrum of bound states of heavy QQ systems. 
The level shift depends on the parameter Ar, where r = A/ttiqu 2 is the typical time of the low lying levels of the 
system, and A is the correlation length of the gluon correlator A defined as D(x — > oo) ~ (: F 2 :) exp(— |x|/A). Thus, 
at large distances the physically motivated asymptotics of the correlator is exponentially decreasing. The gluon 
correlators are the base elements of the stochastic model of vacuum || and in the description of high-energy hadron 
scattering Q. 

Measuring the correlation length and vacuum field correlators was the motivation to investigate these quantities 
on the lattice. New high-statistical LQCD measurements of the gauge - invariant bilocal correlator of the gluon field 
strengths have become available down to a distance of 0.1 fm 0. Recently, the field strength correlators have also been 
studied in the effective dual Abelian Higgs model in || and QCD sum rule approach |J. In all these approaches (see 
also fief] ), the exponential decay of the correlators at large distances has been observed. However, these investigations 
still omit a small and intermediate distances behavior of the nonlocal condensates. 

On the other hand, in QCD there is an instanton flffl , a well known nontrivial nonlocal vacuum solution of the 
classical Euclidean Yang - Mills equations with finite action and size p. The importance of instantons for QCD is that 
it is believed that an interac ting in stanton ensemble provides a realistic microscopic picture of the QCD vacuum in 
the form of instanton liquid [O 13 1 (see, e.g., a review [Q). It has been argued on phenomenological grounds that 



2 



the distribution of instantons over sizes is peaked at a typical value p c ~ 1.7 GeV^ 1 and the liquid is dilute in the 
sense that mean separation between instantons is much larger than the average instanton size. 

In our previous work p5[ , we have shown that the instanton model of the QCD vacuum provides a way to construct 
nonlocal vacuum condensates. Within the effective single instanton (SI) approximation we have obtained the expres- 
sions for the nonlocal gluon Dj"' p<7 (x) and quark Mj(x) condensates and derived the average virtualities of quarks X' q 
and gluons X 2 in the QCD vacuum. It has been found that due to specific properties of the SI approximation (self- 
duality of the field) the D\ gluon form factor is exactly zero. The behavior of the correlation functions demonstrates 
that in the SI approximation the model of nonlocal condensates can well reproduce the behavior of the quark and 
gluon correlators at short distances. Really, the quark and gluon average virtualities, defined via the first derivatives 
of the nonlocal condensates Mj(x 2 ), Dj(x) at the origin, 



2 _ 8 dM^x 2 ) 
q = ~M,(0) dx 2 



p 2 ' 9 dx 2 

are connected with VEVs that parameterize the QCD SR, 



24 1 



X= ° 5 pi 



(3) 



A «- (:qq:) ' ^~ (: F 2 :) ~ 2 (: F 2 :) (: F 2 :) " W 

Where (: fF 3 :) = (: f abc F^F* p F^ :), J 2 = J° J£ and J£ = q(x)^f,q(x). The value of X 2 q » 0.5 GeV 2 estimated 
in the QCD SR analysis 0Jl^] is reproduced at p c s» 2 GeV -1 . This number is close to the estimate from the 
phenomenology of the QCD vacuum in the instanton liquid model. The model provides parameterless prediction 
for the ratio X 2 /X 2 — 12/5 and so, X 2 ~ 1.2 GeV 2 . In (lj| the effect of the inclusion of Schwinger exponent into 
semiclassical calculations has been analyzed. For some quantities this effect is very strong, being of an order of 100% 
for gluon and quark average virtualities. 

Nevertheless, the SI approximation used evidently fails in the description of physically argued distributions at large 
distances. In asymptotics, we have found Mi(x — > oo) ~ pt/x 2 and Dj(x — > oo) ~ Pc/ 2 - 4 f° r the quark and gluon 
correlators, respectively. Thus, the SI solution over mathematical vacuum provides too slow asymptotics at large 
distances. We should conclude that in order to have a realistic model of vacuum correlators, the important effects of 
instanton interaction with the long - wave vacuum configurations have to be included. 

The key point in the picture of realistic instanton vacuum is the interaction of pseudoparticles in the vacuum. 
In [ p^| , the interaction of a SI with an arbitrary weak external field has been examined and dipole-dipole forces in 
a far separated instanton-anti-instanton system derived. Later in j|(J, this background field has been interpreted 
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as a field of large-scale QCD vacuum fluctuations, and the influence of the quark and gluon condensates on the 
instanton density has been considered. The main assumption of the instanton liquid models |B] is the dominance 
of the instanton component in the vacuum and that, in particularly, the gluon condensate is saturated by weakly 
interacting instanton liquid. In deriving instanton ensemble properties the instanton-anti-instanton interaction at 
intermediate separation start to play the crucial role in stabilization of the liquid ]l3| ] . However, it turns out that the 
final result strongly depends on the field ansatz for an instanton-anti-instanton configuration |2l]Jl^] . Further, in all 
instanton-anti-instanton ansatze suggested the influence of the physical vacuum on an instanton profile function has 
not been taken into account and the profile has only power decreasing asymptotics, which contradicts the expectations 
concerning the vacuum field correlators. Moreover, it is known that the instanton liquid is not responsible for large- 
scale effects like confinement |l4j| . Another point is that the instanton density n c in the instanton liquid models is 
normalized by the value of the gluon condensate ^: 2ja F^ V F^ V :) = 0.012 GeV 4 obtained in |lj from an analysis of 
charmonium spectrum. More recently in p2[ , a detailed analysis based on charmonium, bottonium and heavy-light 
mesons have led to a twice larger value of the gluon condensate (: ^ F^F^ :} = 0.023 GeV 4 . Indcfiniteness in 
the normalization provides a window for existence of a large-scale field component in the QCD vacuum along with 
short-scale instantons. 

In the present work, assuming dominance of the weak interacting instanton liquid in the QCD vacuum, we suggest 
that there is also a weak residual component of the vacuum field with a large correlation length R of the order of 
the confinement size. We are going to show, assuming only very general properties of a weak large-scale vacuum 
field, that it deforms an instanton at large distances leading to exponentially decreasing asymptotics of the instanton 
profile and the instanton induced vacuum field correlators. The vacuum model considered is a two-phase one. The 
large-scale phase is described by the background field and dominates at distances compared with the confinement size. 
The short-scale phase is dominated by instantons and is responsible for the spontaneous breaking of chiral symmetry 
and the solution of the U a (1) problem. The vacuum model suggested reveals a definite non-locality mechanism 
in the framework of QCD. We shall illustrate this by analyzing the vacuum gluon form factors D(x) and Di(x). 
Unfortunately, the normalization of contributions from different phases to the gluon condensate is not fixed by the 
instanton model and remains as a free parameter, but the form of the correlation functions can be described in detail. 
The latter is the main motivation of the present work. The determination of the vacuum field correlators is one of the 
main tasks of the theory in describing the non-perturbative dynamics at large distances compatible with the typical 
hadron size. 

The paper is organized as follows: in Section 2 the instanton field in the background of weak large-scale vacuum 
fluctuations is considered; solutions of the vacuum field equations for small and large distances are analyzed separately 
and the constrained instanton solution interpolating two asymptotics is suggested in an ansatz form. In Section 3, the 
space coordinate behavior of the nonlocal correlator of the gluon field strengths is found and the main asymptotics 
of the correlators D(x), D\(x) at large distances are derived. These asymptotics are driven by the strength and 
correlation length of the large-scale vacuum fluctuations, rather than by their form. 



II. CONSTRAINED INSTANTONS IN QCD VACUUM 

The classical Yang - Mills equations in the Euclidean space 

DpF» v {x) = 0, (5) 

where the covariant derivative is = — iA^Ta/2, have an instanton (+) ((— ) is for an anti - instanton) as a (anti) 
self-dual finite-action solution with topological charge, Q = ±1: 

2 

A-ting,^ x o) = lOf-rf^ (x - x ) v </ g {x - x ) with <p g (x) = ^ (in singular gauge), (6) 

X yX — (— p J 

localized in size p. In (^|), xq is the position of the instanton center, Oj b is the ortogonal matrix of instanton 
global orientation in the color space and = tAa^u T —^abc^-bc^v ar< 3 t'Hooft symbols. The solution (||) is written 

in the singular gauge within the SU{2) subgroup (with generators I f ) of the SU C (3) theory. This classical field 
configuration reflects the symmetries of the initial theory in terms of collective variables corresponding to translational 
transformations, rotational symmetry in color space and conform transformations. 

The solution (^) is given in mathematical vacuum and has an unpleasant property of a very slow decay at large 
distances noted in the introduction. This situation is inadequate since the physical vacuum is not empty but looks like 
a medium densely populated by large-scale vacuum field fluctuations. In the background of the large scale fluctuations 



3 



there are developed non - perturbative fluctuations of a smaller size among those instantons which play a dominating 
role. The long-wave gluon vacuum field, which is the background for a selected instanton, may be of a more general 
origin and phenomenologically can be parameterized by the vacuum correlation functions of the gluon operators 
contributing to the corresponding nonlocal condensates. 

What is important is that at random (stochastic) background vacuum field with fixed vacuum expectation values 
of singlet operators the scale invariance of the effective theory is spoiled already at the quasi-classical level and the 
instantons are no longer exact solutions of the field equations and the Dirac operator has no zero modes. 

The similar situation has been observed in the standard electroweak Yang-Mills-Higgs model. There, the background 
of Higgs field with nonzero vacuum value (tp) and coupling A affects the instanton configuration J23|] . In the presence of 
(even small) effects violating scale invariance of the initial theory, instanton solution does not exist at all. Nevertheless, 
as it was stated in |2j| and fairly explained in pi|-p6[ , if the Higgs field is rather weak and some additional constraints 
are introduced, there may be constructed an approximate solution, so-called, constrained instanton (CI). These 
constraints limit the degree of freedom along the size p parameter. The constrained solution at small distances 
|a;| << VA (ip) 1 approximately has the form of an instanton, and at large distances |a:| >> VX (p)^ 1 >> p has the 
asymptotics of massive (gauge boson) particle exp(—g (tp) \x\/\/X). In p4| , it has also been noted that the gauge field 
propagator of the CIs decays exponentially at large |a;| and thus does not affect the long-range behavior of the theory. 

The aim of this section is to show that an analogous phenomenon takes place in QCD considering an instanton field 
A^x) in the physical vacuum. In distinction with the Yang-Mills-Higgs model, in pure QCD there is no Higgs field 
from the beginning and it is the long-scale vacuum gluon field, b^(x), that models a source perturbing an instanton 
at large distances. The deep reason of this effect is in an existence of the quantum anomaly in the trace of energy- 
momentum tensor |p7| . It will be shown that the presence of this background field, characterized by its vacuum 

expectation value (^(f^^ y and correlation length i? (introduced below), sets the final scale and defines the 

r// s2\ 1-V3 

>> p. Here and below, 



b 

-1/3 

deformation of the instanton solution in the asymptotic region \x\ > 



nv) > R 



b 

(■■■)& = Jdo~[b] ... means the average over nonperturbative random background field weighted with some measure 
da [b] . This solution is stable against shrinking the instanton to a point if some constraints are added. By analogy 
with the solutions analyzed in the Yang-Mills-Higgs model |24j , we shell call these interpolating fields constrained (or 
deformed) instanton solutions. 

In the following we analyze the vacuum field configuration of the single (constrained) instanton A^ 1 (x) of fixed 
size and orientation in the color space in the background of the large-scale topologically neutral random vacuum field 

A„{x) = A^ I (x,x )+b l ,(x), (7) 

with gauge transformation property 

A^x) -> CTt (x) A^x)U (x) + ^[/ t (x) d^U (x) , (8) 

where U (x) is a gauge transformation matrix. Considering the field ansatz (|7|), one has to take the instanton field 
A^ 1 (x) in the singular gauge ||l3| , |l4"f 

A C £t;}{x) = 20f V b >± (x - *„)„ p c g \x - x ), (9) 

*W*)U, = 1 ' ^)l™ = °- 

The last conditions mean that the constrained instanton has a finite action and a modulo unit topological charge, 
but, in general, ceases to be self-dual field. In the coordinate space the instantons in this specific gauge fall off rapidly 
enough to provide a weak interaction with the background field and the quasi-classical approach is justified. The 
weakness of the interaction allow us in the following to neglect the back reaction of the instanton on the background 
field. As to nonperturbative background field it is convenient to choose it in the Fock-Schwinger gauge [ p8| 

(x - xo)^ b» {x - x ) = 0. (10) 

In the following we put the instanton center Xq to the origin of coordinates x = 0. 



2 The similar model has been considered earlier in [ ^7| . 
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As an illustrative model for the background field, one can keep in mind the self-dual homogeneous vacuum gluon 
field b a ^ (x) = \n a b llv x ll , b^b^y — b 2 , where n a is the orientation vector in color space and b^ v is the constant 

field-strength tensor |30|. The corresponding measure / da [b] — / °° dbD (b) J d£l f dft c averages over field amplitude 
and its orientations in configuration and color spaces. This field configuration with infinite correlation length R = oo 
and infinite topological charge quite correctly describes situation at small and intermediate distances comparable with 
instanton size, but at larger distances the effect of finite correlation length of physical background becomes important. 
Introduction of finite correlation length can be imagined as the inclusion of domain structure in the vacuum p9| . This 
kind of considerations are in the base of the stochastic vacuum model 0. In the absence of a consistent theoretical 
approach to the large distance dynamics one is led to elaborate the problem phenomenologically. 

We assume that at small distances the CI field dominates and the background field b^ (x) is regarded as a pertur- 
bation on A^ 1 (x). At distances much larger compared to the instanton size p the background field b^ (x) is still weak, 
but strong enough to deform and suppress the instanton field. 

The field strength can be written as 

f;„ [a ci + b} = F™> a + f & v + af;„ [a ci , b] , (11) 

where = F^[A CI ], F^ v = F£ v [b] and 

AF;„ [A CI ,b] = f abc (A c /Hl + A° v ^b% 
and the effective Euclidean action of the instanton in the random background field becomes 

{F™>*F°!" + AF£ V [A CI , b] AF; V [A ci , b]}^ (12) 

In deriving these expression we have used the color-singlet properties of the large-scale vacuum on an average 

(*M„>6 = ( 13 ) 

and average over relative orientation of instanton and background field in color space. We also neglected the terms 
higher order in interaction. Below, these terms will effectively be accumulated in the form of constraints below and, 
what is important in the present consideration, they do not influence the form of the solution at asymptotics. 

Similar to Affleck analysis p3|, we come to the conclusion that for a background field configuration b fJi (x) with a 
fixed nonzero condensate value no instanton solution exists. This can be seen from the rescaling x — > ax, A*? 1 (x) 



— > a 1 A^ I (ax), b )1 (x) — ► b^(ax), preserving finite vacuum average (F^^j — const, under which Se transforms to 

s E - ^ (J d 4 x {f^fc 1 '* + «~ 2 af; [a ci , b] af; v [a ci , b]}^ (u) 

If A^ 1 (x) is a stationary point, then dSs/da = and the action is minimized by the instanton of vanishing size. 
Thus, given any field configuration we can always rescale it to get smaller action, except in the trivial case. 

Now, let us consider the problem from the point of view of the equations of motion for the deformed instanton in 
the background of large-scale random vacuum fluctuations, that follows from (|l2|), 

Df [A CI ] F c J- h + f bac f bkl (A c ; j ' k (b%) b - A c v l - k (b%) b ) = 0. (15) 

In the Fock-Schwinger gauge the background field has a representation in terms of its strength 

b^(x) = [ daaF£^{ax)x v (16) 
Jo 

and the bilinear field averages become 

(b«(x)b b v (x)) = / da [ d(3a(3x p x a (F b a ^(ax)F b b ^(l3x)).. (17) 



o Jo 



In the non-Abelian case the correlator in the integrand of (|17j) is not gauge- invariant, however, in the Fock-Schwinger 
gauge this correlator coincides with the gauge-invariant correlator in which field-strengths are connected by the 
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Schwinger phase factors E(ax, 0)E(0, fix) in the adjoint representation. Thus in this specific gauge the gauge variant 
left side of (0) can be expressed in terms of gauge invariant quantity ^F£ pfl (ax)E(ax, l3x))F bjJU {flx)^ . Due to the 
gauge-invariance the latter correlator admits a physically motivated model 



x p x a {F a m {ax)F b av {f3x)) = 



=x(a-/3) 



where the function 



B (z 2 ) = D (z 2 ) + D 1 (z 2 ) + z 2 dD 1 (z 2 ) /Oz 2 



(18) 



(19) 



is defined via the form factors D (z 2 ) and D\ (z 2 ) parameterizing the gauge-invariant two-point correlator (||) of the 

background field strengths, with normalization D(0) — k, D\ (0) = 1 — k. The contribution of the background field to 
the gluon condensate is denoted by {F 2 ) b - With these definitions the equations of motion of the CI field interacting 
with a random large-scale vacuum fluctuation field, (|l5|), can be cast in the form 



d; 1 



[A CI ] F^ (x) 



24(7V 2 - 1) 



x 2 $ (x 2 ) Af/> a (x) = 0, 



(20) 



where 



$ (x 2 ) =4 da d(ia(3B {a - (3) x 



$(0) = 1, 



(21) 



and iV c is the number of colors. 

Let us discuss the properties of the solution of Eq. (p0|). In the absence of the background field {F 2 ) b = there 
exists an instanton solution (^|). For (F 2 ^ small enough, such that (F 2 ^ <C l/p 4 , we should expect to find a solution 
of (|2^) in perturbation theory in small parameter (F 2 ) b p 4 , which reduces to (||) when {F 2 ) b — > 0. However, such a 
pcrturbative solution does not exist. The reason is that for the higher order in perturbation terms appropriate finite 
action boundary conditions at large distances cannot be enforced |24j ]. 

The operators that act on higher order terms possess zero mode dAz 1 /dp 



OA 



ci 



dp 



dp J 



r)A CI 



(22) 



which determines a priori the behavior of perturbative around instanton terms at infinity. A way of getting around 
this difficulty J24J] is to extremize the action Se, subject to a constraint. The choice of an explicit form of the 

constraint is quite arbitrary. In p3 it has been proposed the global constraint of the general form 



(A)= d i x[O d (A)-O d (A CI )] =0, 



C7i 



where the gauge- invariant local operator Od(A) has a canonical dimension d > 4. The relevant stationary configuration 
will be a solution of the equations of motion d2G) but with the constraint term added into the right-hand side 



SS E {A) 



SA„ (x) 



5C cons tr {A) 



A CI 



SAv (x) 



A CI 



(23) 



The Lagrange multiplier a in ( p3| ) is to be determined order by order in perturbative theory in (Fjf^p 4 , which 
provides the correct boundary conditions for the higher order terms. The constrained instanton A is the unique 
solution of ( p3|) obtained by this procedure, the j4 c/ -solution turns to (||) when (F 2 ) b p 4 — ► 0. Unfortunately, this 
kind of constraints is nonlinear and the higher order terms in (F 2 ) fe p 4 , depending on the constraint, are difficult to 
evaluate in practice. 

Another way has been suggested in p6[ , where a linear constraint of the general form 



CI 



(A) = / d 4 xtr { (A, (x) - Af (x)) ft (x) } = 



(24) 



G 



has been introduced. It has been proposed, p^j , instead of fixing the constraint to solve the equation for A? 1 (x), 
which is almost an impossible task, to choose A^ 1 (x) first and then find the constraint fP (x) itself by substituting 
A^ 1 (x) into the left hand side of the constrained equation (^3j) . In this way, the freedom in choosing the constraint 
can be used to find it by a given solution. 

One can also restrict by considering only the local operators defining constraints that fall down at infinity more 
rapidly than the interaction term in Eq. (^Tj|) . Under this condition, it is easy to obtain the behavior of the instanton 
enveloped in the background field at distances far from the instanton center. This large-distance asymptotics of CI, 
like its behavior at small distances A 1 ^ (x; p), is a constraint-independent part of the solution. 

We are interested in the asymptotic behavior of the function $ (x 2 ), Eq. (pl|), where the interaction term becomes 
dominant over the instanton self-interaction. It is nice that the leading asymptotics of $ (x 2 ), $ (x 2 ) ~ R/\x\, where 
R is a correlation length, is independent of a particular form of the function B(z 2 \ This property is due to a specific 
dependence of the argument z 2 = (a — (3) 2 x 2 . Let us illustrate this property using a few natural ansatze for this 
function 



B M (x 2 ) = R 2 / (x 2 + R 2 ) , B E (x 2 ) = exp(- \x\ /R), 
(Monopole) (Exponential) 

that provide, respectively, (see appendix B) 



B G (x 2 ) = exp(-x 2 /R 2 ), 
(Gaussian) 



(25) 



r s (x 2 ) 



2\ \x\— >oo 8 R 



Q a *T" I' 
3 \x\ 



where 



a$ = 0(1) = 



1 



monopole, 
exponential, 
gaussian. 



(26) 



The exponential ansatz (modulo powers) appears as an asymptotics of the solution of the Yang-Mills-Higgs model 
p4] , p5 | and also is used in parameterization of large-scale behavior of Lattice QCD data [Q. The monopole form 
resembles the asymptotic behavior of the SI correlator pq| , where, in fact, it has steeper decrease like 1/x 4 . 

Asymptotic behavior of the instanton solution deformed by large-scale vacuum fluctuations at large Euclidean 
distances |ar| — > oo can be derived from the analysis of the equation 



. 1 

a<$N c 



0. 



Vg = 



9(N 2 



which follows from Eq. 20. Due to a decreasing character of the field asymptotics, A^ 1 — > at 



(27) 
(28) 

oo, only linear 



terms in the short-wave CI field A^ 1 are kept in the kinetic part of equation 
defined by (77JJ = 77^° in the following) 



For the profile function (p^ s (x 2 ) 



ACI,a( \_—a ZJL ln as 



(29) 



we find from the asymptotic equation (p7|) the large distance solution 



*7 (x 2 )~K 



4/3 



(Vg \x\) 



.3/2 



(30) 



where K v {z) is the modified Bessel function with index v having the asymptotic behavior K„(z) — * \/j^e~ z as z — > 00. 
We have to note that in the case of the homogeneous background field with infinite correlation length one get the 
equation similar to ( B7I), b ut with the coefficient proportional to x 2 in the last term that results in the Gaussian form 
of asymptotics ( c.f. , ( |27|] ) ) . 



3 At this point we have to note that the influence of the instanton ensemble on the instanton profile has been discussed in fL3| , 
The authors found that this interaction perturbs the self-interaction part by the term /j^A^. It turns out that numerically the 
coefficient p. 2 g strongly depends on the instanton-anti-instanton ansatz chosen J2TJ and, as it is seen from (|2^), has subleading 
behavior in the limit of large x. 
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The requirement of the instanton-background interaction being weak allows us to consider the dimensionless pa- 
rameter a g = rjgp a small value. It means that the region, where the instanton field dominates (small distances), and 
the region, where the background field dominates (large distances), are well separated and the large distance effects 
do not destroy the instanton. Then, the overall constant is determined by matching, at distances p <C \x\ <C Vg^ 1 , the 
leading terms of the expansions of A p 2 (x) at small distances, which is an instanton A 1 (x), and at large distances, 
which is an asymptotics (p|), (J30|) , 



A c /^(x)=rf^a 4/3 a 2 g K i/3 



(Vg \x\) 



3/2 



where 



«4/3 



r(l/3) 3V3 



(31) 



is the normalization coefficient and T(z) is the Gamma-function. 

Thus, we find that in the singular gauge the CI solution has an exponential decreasing character (|3l| ) far from the 
instanton center. This behavior sharply differs from the power decreasing asymptotics of SI (^|). This modification 
of behavior follows from the fact that the instanton solution is considered in the physical vacuum populated by the 
large-scale gluon field fluctuations. The background field modifies the long-distance behavior of the instanton and 
leads to appearance of the "second scale" parameter A g = l/rj g (see |l5[) in gluon distributions. At the same time, 
the effect of long- wave vacuum fluctuations is not very essential for the behavior of the instanton at short distances. 

Now we are looking for the constrained solution of (GCl) in the ansatz form 



A? a (x) = 2rf^x vV>g (x 2 ) 



which has the behavior at small and large distances 



(32) 



(x 2 ) = 



(x 2 4 

a 4/3 a l K 4/3 



+ (%M) 3 y>(i) (rr 



(Vg \x\) 



3/2 



P 



at 



<pW (r, g \x\). 



M-o, 

at \x\ — > oo, 



(33) 



where the forms of the expansions can be deduced from expanding the leading terms with respect to j^y and r\ g |a;|, 
respectively. The first terms in the expansions are exact and constraint-independent ones, however, as it was shown 
in 24|], all higher order terms have dependence on the constraint chosen. 

The condition of finiteness of the action, which is also constraint-independent one, should be imposed on the 
desirable solution. To this goal, let us rewrite the CI field strength as 



F% Ia (x) = 4 [%„ux (x) + (x^ p - x v % p ) x p u; 2 (x)] 



(34) 



with forms 



wx (x) = x 2 if 2 g (x 2 ) - (fig (x 2 ) , uj 2 (x) = f 2 g (x 2 ) + 



difg {x 2 ) 



to be used further. With this parameterization we have the action 

1 



dx 2 



d^x [f; v {x)f; v (x)] 



CI 



(35) 



(36) 



with the action density 



(x) F« (x)]] CI = 96 K (a;) + u; 2 , (x)] 



(37) 



where 



LU 3 (x) — X 2 UJ2 (x) — U>x (x) ■ 

Now, if we use the singular gauge for Ar^ 1 , then, to guarantee finiteness of the action, the condition has to be fulfilled 



V fl (x 2 ) |x^o ^1 + (x 2 ) . 



(38) 



For further references we present here the well-known expressions for the SI profiles in the singular and regular gauges 
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< l9J (* 2 ) = 4s, W (* 2 ) = i (39) 



x x p 



sing,I ^ = £_ x 2 u s 2 m3j ( x ) = 2wf" 9 '- r (x) , (40) 

( x p) 

ol^( x )=-f ^• I (x) = 0, (41) 

[ x p) 



and gauge-independent expressions for the density action and the action itself 



[^(*W*)]'=7^£, S£ (42) 

By using the asymptotic properties of the CI solution ([33]) and finite-action condition (pq) , which are constraint- 
independent, we are able to construct ansatz. Certainly, this procedure is not unique and in principle one can impose 
further physical requirements to constrain the behavior of the solution in the intermediate region. These details, 
however, can be taken into account by choosing proper constraints. Thus, the freedom in choosing the constraint can 
be used to find it by a given solution, instead of solving complicated equations p(| . 
Let us consider the following ansatze for the CI profile written in the singular gauge 

, 2 , 1 K4/3 [z PtX ] 



4/3 [Zp,0 



where we have introduced the notation 



(* 2 ) = U^fX 2* > (4 5 ) 
x z [x z + p (x z )) 

^ 3/2 / 2 1 2\ 3 / 4 
Z P-X = ^Vg [X + P ) , 

p 2 (x 2 ) = a 4/3 a 2 g x 2 K 4/3 [z , x ] , p 2 (0) = p 2 . 

Note that all ansatze have easily identifiable instanton parameters. By translational invariance the center of CI can be 
shifted in ( (43| ) - ( p5| ) from the origin to an arbitrary position xq : x — > x — xq . The CI profile functions corresponding 
to these ansatze are shown in Figs. (|l|) and @ along with the instanton profile (g9|). Figs. (|l|) and (||) display 
the dependence of the profile function on the external field parameter prj g and on the form of ansatz, respectively. 
To make the difference more clear, we also take for illustration a large value of the parameter {prig) 2 = 3. We see 
that if at small distances the CI is close to the instanton form, then at large distances this solution has exponential 
asymptotics instead of power-like for the instanton. 

Last two ansatze are similar to ones suggested in p6| , where the preference has been given to the (p 3 (x 2 ) form, 
since it has better convergent properties in the expansion of CI (|33]). Moreover, with this profile the constrained 
solution in the regular gauge looks similar to the SI case 

^ CI (x 2 ) = - ' 



+ p 2 (x 2 ) 

In order to pass from the constraint instanton in a singular gauge to the instanton in a regular gauge one can translate 
the general gauge transformation (||) into the form 

A^{x) -> & (x) A° x (x)Sl (x) + iflt (x) d„n (x) , (46) 



b^x) fit ( x ) b„(x)Sl(x), 

with transformation matrix 







We have numerically calculated the dependence of the CI classical action, Eqs. ([36|) and ([37]), on the instanton size 
p. This dependence for the three ansatze (0), and ( ^5|) is shown in Figs, (g) and (||). We have to stress, that the 
profiles of the field A^ 1 and the action S^^iepcnd on the choice of constraint. However, the full effective action, with 
the terms coming from the Jacobian included, is constraint-independent [2(J. The weak dependence of the action, 
Sfi 1 , on the profiles ^1,2,3 (x 2 ) (see Fig. (^J)) indicates that in the region of parameters pj] g < 1, the influence of these 
additional terms is small and the exponential part of the action, S^, can be used as a good approximation. We see 
in Fig. (^J) that the CI action is larger than the instanton one and monotonically grows with the instanton size. It is 
natural because the CI- "solution" does not self-dual one and does not realize the minimum of the action. Instead, it 
represents the bottom of the valley parameterized by the quasi-zero mode p. 

We are not going to discuss further details in construction of the total effective CI action, that takes into account 
small quantum oscillations around the nonperturbative configuration (Q) and the interaction of constrained instantons, 
and postpone them until further publication. Just point out that other effects dominating the effective action at small 
p come from the running coupling constant g 2 (p) and the path integral measure over the size of instanton dp/p 5 . It is 
well known that in the model of the coupling constant, which freezes it to a constant at some large po, the corrected 
action 

Sl OT (P) = J d*x [F% (x) F- v (Z)] CI + 5 In p, 

as a function of the instanton size, has a minimum. The position of the minimum is correlated with the freezing 
parameter pg, which can be chosen to provide the value p m in ~ 2 GeV -1 , and the environment of the large-scale 
vacuum fluctuations makes the minimum more prominent (see ifjlf for recent discussion of similar results). However, 
the typical CI action at p m in is rather large, its numerical value being around 25. This means that the configurations 
with small number of instantons and anti-instantons are not important statistically and suggests that the interacting 
instanton and anti-instanton ensemble could be a more important type of configurations. The leading interaction 
term of a widely separated instanton - anti-instanton pair in physical vacuum as described by Eqs. ( f43| ) - J45| ) 

falls with separation L like exp A/3(r) g L) 3 ^ 2 ^j and differs from power-like decreasing behavior found in (lj]jL| in 

unconstrained case. In the following, considering the nonlocal properties of gluon condensate, we accept that the 
instanton liquid is formed due to the instanton-anti-instanton interaction and the instanton density n c and size p c are 
fixed @:n c wl fm" 4 , p c ^ 1/3 fm. 



III. SHORT-RANGE VACUUM CORRELATORS IN THE CONSTRAINED INSTANTON MODEL. 



Within the model considered the full gluon correlator may be conventionally written by using the expression for 
the field strength (O) as 



<: F^[A CI + b] (x) F P M CI + b] (y) :) = {F% (x) F C J (j,)) + 



+ 



F btllv (x) F b , pa (y) :) + (: AF^[A CI , b] (x) AF pa [A CI , b] (y) :) mterf , (47) 

where the brackets ( ) mean averaging over vacuum fluctuations (^) and we do not display Schwinger phase factors 
explicitly. The last term represents the interference of short- and large-scale fields and will be discussed below. For 
the large-scale correlator (|l^) we have already suggested the model for the form factor B{x 2 ) in (p5|). Now, let us 
calculate the short-range part of the gluon correlator. 

Let us construct the correlator D^' pcT (x — y) of gluonic strengths (|l|) in the quasi-classical approximation by using 
the CI solutions given by Eqs. (|32] ) and (^3|) - (|4g). We will use a reference frame where the instanton sits at the 
origin and a relative coordinate (x — y) M with respect to the position of the instanton center is parallel to one of the 

coordinate axes, say fi = 4, serving as a "time" direction (i.e.,x— y= 0, 2:4 — 2/4 = \x — y\), and reduce the path 
ordered exponential to an ordinary exponential 



E(x, y) = Pexp 



A^(z)dz^)=LHx)L(y) 



(48) 



with 



L (x) = exp =Fz t 



-«( 



,X 4 



= T«r± • x» (x) . 



(49) 
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where 















X 


,x 4 ) = 


X 


Jo dtif>9 ( 


X 





(^i,f), x° (x) = cos a (x) , x J (x) = ^x 1 / x^sina(x) 



(50) 
(51) 



The factor L (x) coming from the Schwinger exponent can be accumulated in the definition of the field. This repre- 
sentation of the field may be called the axial gauge representation A^zjn* 1 = 0, since in this gauge with the vector 



Tin Xt, 



the Schwinger factor E(x,y) = 1. 



In the CI background the bilocal gluon correlator acquires the form 
D» v >P°(x) = <: Tr (F {ax) ^(0)F {ax)pa (x)) :>= 



d 4 z 



dQTr ( F? , (z - -)Ff , (z + -) 



(52) 



where nf is the effective instanton / anti - instanton density, z is the collective coordinate of the instanton center and 
fl is its color space orientation. 

To extract form factors D (x 2 ) and D\ (x 2 ) it is easier first to average over the instanton orientations in the color 
space and take the trace over color matrices by using the relations 



Nr. 



d u bi 



(53) 



T a T b = S ab + i£ abc T c . 



V 'v — ' 0, lfiU 1 ' 

Then, it is convenient to define the combinations of functions D (x 2 ) and D\ (x 2 ) 
A (x 2 ) = S^Su,, \' OI = D (x 2 ) + £>i (x 2 ) + -x 2 



B (x 2 ) = 4 ^<W-^-^4t = £> U 2 ) 



<°l^l°) 



<9x 2 

. 9 clDi (x 2 ) 



(54) 



taking the boundary condition, D(Q) + -Di(O) = 1 and the asymptotic conditions D(oo) = Di(oo) = 0. After direct, 
but cumbersome calculations we come to the expressions for the functions A and B: 



A(x 2 ) = -N 



D 



drr 



dt{[u>i (z+)u>i (z-) + uj 3 (z + )u> 3 (z-)] (3 -4 sin 2 (a*)) - 



(55) 



-2uj 2 (z+) UJ2 (z-) [r 2 x 2 (1-2 sin 2 (a z )) - rx (z+ ■ Z-) sa\(2a z )] } 



-1 r* pOC pOO 

B(x 2 ) = —N D drr 2 dt {ui 1 (z + )uji (z_) (3 — 4sin 2 (a z )) - (56) 
n Jo Jo 

-uji {z+)lu 2 (z-) [z 2 _ + 2t 2 _ (l - 2sin 2 (a z )) + 2ri_ sin(2a z )] - 

-ui 2 (z+) uji (z-) [z\ + 2t\ (1-2 sin 2 (a z ))) - 2rt + sin(2a z )] + 

+uj 2 {z+)uj 2 (z-) [z\z 2 _ + 2t+t_ (z + ■ z-) (1 - 2sin 2 (a z )) + 2rxt+t- sin(2a 2 )] } , 

where z± — (r,t±), t± = t ± — , the forms lu\ (z) , u>2 (z) , and ^3 (z) are defined in (j35|), No is the normalization 
factor 



N^=G dyy 3 (u> 2 (y) + u 2 (y)) 
Jo 



(57) 



and the phase factor 



2 drcpg (r 2 + (t + r) 2 ) 



11 



reflects the presence of the E exponent in the definition of the bilocal correlator. These expressions for the field 
strength correlators are general for any field given in the form (|32|). The gauge invariance of the functions A{x 2 ) 
and B{x 2 ) can explicitly be checked for Exp.(|55|, Kq) by transforming, for example, the field from the singular to 
regular gauges, Eq. (|4g). The expressions for A{x I ) and B(x 2 ) may be considered as generation functions to obtain 
condensates of higher dimensions in the instanton model approach. From a technical point of view this procedure is 
more convenient than their direct calculations [p2|]l5| |. 

In the SI approximation the form factors BUx 2 ) — A 1 (x 2 ) reproduce the expression, Eq. (21) from |l5), for the 
gauge-invariant correlator. As it has been shown in |l5f| (see also p4j| ), in the SI approximation the term with the 
second Lorentz structure, Di(x 2 ), parameterizing the gluon correlator (^ ) does not appear. This fact is due to the 
specific topological structure (self-duality) of the instanton solution. Both the Lorentz structures arise in the r.h.s. of 
(|52|) if one takes into account the background fields. 

The form factors D{x 2 ) and D\{x 2 ) are determined numerically by solving the equations ([54]) and plotted in Fig. 
(|J) in coordinate space and in Fig. (|6|) in momentum space. The constant k defining the relative weight of D 
functions (D (0) = k, D\ (0) = 1 — k) depends on the background field strength parameter rj g p and it is close to one 
in the region of reasonable physical parameters. It is equal to k = 1 at (rj g p) 2 = (SI case), k = 0.997 at {r/ g p) 2 = 0.1 
and k = 0.926 at (r/gp)" = 3. In all the cases the parameter n is close to one in accordance to the fits of lattice data 
p3| . It means that in the weak background field the dominant role of the D{x 2 ) function remains as in the SI case; 
it is close to the SI form at small distances and exponentially rapidly decays at large distances. The Dx(x 2 ) function 
is small and positive everywhere, its behavior is very sensitive not only to external field but also to the gauge phase 
factor effect^. In Appendix A we show that the reason for smallncss of the Di{x 2 ) function is "almost" self-duality 
property of the CI solutions. Both the form factors possess two zeros at large distances and at very large x develop 
positive asymptotics 

Dix 2 )^^ 2 )- |zr 3 / 4 ( % p) 4 exp (-0.473 (% M) 3/2 ) , (58) 

where r\ g ^ and the constant in the exponent is found numerically 

In order to have contact with the QCD vacuum phenomenology and specify further the instanton-induced model 
of the gluon correlator, let us discuss the contributions of different terms in p7| ) to the gluon condensate 

(: F^[A CI + b] (0) F, V [A CI + b] (0) :) = {F™ (0) F™ (0)) + (59) 
+ (: F btia , (0) F M „ (0) :) + (: AF, V [A CI + b] (0) AF^[A CI + b] (0) :} Merf . 
The background contribution to the gluon condensate 

(:F b ^(0)F b ^(0):) = (F 2 ) b (60) 
serves as a parameter of the model and, by assumption, is much smaller than the CI contribution given by 

(0) tfL* (0)) = 32n 2 n c N D \ (61) 

where N^ 1 « 1 (see Fig. [I]). The interference term after averaging over relative color orientations and using relations 
(|l3|), (pTp acquires the form 

{:F» V [A CI + b]F, v [A CI + b]-) interf = (62) 



where $(^ 2 ) is defined in (|2l]) (the explicit forms of $(z 2 ) are outlined in Appendix B). 

The interference term depends on two dimensionless parameters a g = p c r\ g and (3 = p c /R and the background 

. 4 _9(JVg-l) 
a$iV c 



field condensate can be parameterized as (F 2 ) b p\ — — - — ^— - — a ^f3- The instanton size p c comes in the last formula 



4 In ref. [Q, in the calculations of the form factors D,D\ based on the instanton-anti-instanton ansatz both the influence of 
the physical vacuum on the instantons and the gluon correlators as an effect of the P— exponential factor on the form factors 
has been ignored. As a result, the negative D\ has been obtained in [[ji]. However, both the facts are important in determining 
the correct norm and forms of the form factors, in particular, in obtaining a small Di 
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with high power and leads to indefiniteness of the factor of order 2 in the relation of the external field condensate 
to the parameters a g and (3. To reduce this uncertainty, we can use the physical information about the vacuum 
properties provided by the QCD SRs and lattice QCD. Indeed, as it has been shown in |lE( ], in the SI case there 
is a relation between the instanton size and the average virtuality of quarks in the vacuum, Eq. (|§J). The value 
of the average quark virtuality has been estimated in the QCD sum rule analysis, A 2 = 0.5 ± 0.05 GeV 2 in pi; 



i 

A? = 0.4 ± 0.1 GeV^ in pi, and from the lattice QCD calculations A 2 = 0.55 ± 0.05 GeV 2 in pf. The relations 



(H) remain good approximation in the CI case if the external field does not strongly deform the instanton. Numerical 



calculations of A 2 , defined in (|3D , lead to estimates 



A 2 =4.8^ K=0), A 2 = 5.7-^ (a 2 g = 1). (63) 

Pc Pc 



We show below that physically motivated background field has the strength parameter a g < 1 and thus the value of 
A 2 increases not more than 20%. 

The relation for A 2 in (||) can be used to get the scale for the background field condensate 



and we accept in the following 



, 9N 9(V 2 -l) , 



A 2 = 0.5 GeV 2 . 



What is the expected range for the parameters a g and /3? By analogy with the instanton liquid vacuum, the parameter 
(5 can be interpreted as a ratio of the instanton size to the inter-instanton distance and it is adjusted as (3 ~ 1/3. 
Then, the estimate of the upper limit for the strength parameter a g follows from the assumption that the contribution 
of the background field to the total gluon condensate (0 \ F 2 \ 0) total ~ 1 GeV 4 |^2| is quite small. This assumption 
reduces the influence of the model dependent part of correlator and leads to the bound a g < 1, or for the dimensional 
parameter, tj 2 = a 2 A 2 /4, rj g < 0.35 GeV. In Fig. (Q) we present the values of the interference term as a two parametric 
plot and see that its contribution to the gluon condensate is small if the short-range and large-scale fluctuations are 
well separated: a g < 1 and /5 <C 1. For completeness in Appendix B we present the small interference contributions 
to the functions A(x 2 ) and B(x 2 ). 

Thus, we construct the model of the gluon correlators. Within this model the functions A(x 2 ) and B(x 2 ) are the sum 
of the short-range instanton induced contributions ( |55| ) and ( |56| ) multiplied by the weig ht factor 32?r 2 n c / (0 \F 2 \ 0) total 
and the long-range contribution (|l|) modeled by Exp. @ with the weight factor (F 2 ) b / (0 |F 2 | 0) totar The 
parameters of the model are the average instanton size p c « 0.3 fm, the effective instanton density n c « 1 fin ~ 4 , the 
strength \F 2 ), < 32tt 2 ti c and the correlation length R « 3p c of the background field. The first two parameters are 
estimated within the instanton liquid models, being in the dilute liquid limit expressed through the vacuum averages: 



p 2 c = 2A- 2 and n c = (2irI/N c ) (\(qq)\ 2 /A 2 ) , 



where the numerical constant / m 0.6 |35|| . The latter relation is a 

consequence of the gap equation [[L3|]. The form of the short-range correlator is defined by p c at small distances and 
by the long-scale parameter rj g at large distances. The form of the long-range correlator at large distances can be 
motivated by the results obtained in the dual effective model of QCD || , where they have exponential decrease (modulo 
powers) similar to exponential ansatz in (p5|). Basing on the results of the dual model and lattice measurements one 
can expect that A(x 2 ) « B{x 2 ) for the long-range part of the correlator. 

The field-strength correlators have been studied on the lattice in [Q, [^0|. There, the following two combinations 
of form factors have been measured: 

D ± (x 2 ) = D (x 2 ) + D 1 (x 2 ) , 



D^x 2 ) ^D(x 2 ) +Dl (x 2 ) + x 2 ^l. 



where D± (a; 2 ) = 2A (a; 2 ) — B (a: 2 ) and D|| (a; 2 ) = B (x 2 ) in terms of the combinations defined in (|54|). The 
lattice measurements of the field strength correlators are also obtained with the straight line path in the Schwinger 
exponent. The direct comparison of the model calculations with the lattice data is a delicate problem, since the used 
parameterization is rather conventional in separating the residual perturbative tail (divergent term ~ a;~ 4 ) from the 
non-perturbative part (pure exponential finite term). As a result, the fits with ad hoc chosen parameterizations are 
very unstable with respect to the extraction of the quantities of physical interest: the correlation lengths, the gluon 
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condensate, etc. |3^]. The reason is that the perturbative part is strongly divergent (as it is seen from lattice data), its 
contribution at small distances would be strongly dependent on the parameterization procedure. On the other hand, 
by construction we calculate the non-perturbative part of the correlators with perturbative contributions subtracted. 
In future, it would be quite desirable to make a new fit to the lattice data using the Eqs. (|55|), (156]), as an input for 
a non-perturbative part of the correlators. 

At the present stage, we restrict ourselves only to a few qualitative remarks. In ref. [j33f , the range of values of some 
physical quantities was discussed which can be fitted from the lattice data according to different parameterizations. 
Namely, the correlation length, the gluon condensate and the normalization of form factors have been analyzed. As 
it has been noted above the normalization of form factors k, consistent with small D% (a; 2 ), is in agreement with the 
instanton model and that the value of the gluon condensate serves as a free model parameter. The values of the gluon 
condensate extracted from the fits to lattice data are very sensitive to the parameterization used, being within the 
interval <: ^ F£ V F£ V :>= (0.005 - 0.03) GeV 4 . In the lattice "full-QCD" fit of an average correlation length l G of 
the gluon strength, defined as Q 

1 f°° 

lG = mv)j dxD ( x2 y ( 64 ) 

is in the range Iq ~ 0.35 — 0.45 fm with lattice quark mass am q = 0.01 and Iq ~ 0.3 — 0.4 fm with am q = 0.02, where a 
is the lattice unit. One can expect, following linear extrapolation, that in the chiral limit am q — ► 0, la ~ 0.4 — 0.5 fm. 
Now let us omit an important but unsolved problem about the difference of lattice and CI renormalization schemes 
and norms, to draw rather a rough comparison of the corresponding results. The instanton model predicts for the 
same quantity: Iq = 0.43 fm at p c r\ g = (SI), Iq = 0.37 fm at {p c Tjg) 2 = 1, — 0.31 fm at (p c Vg) 2 — 3- Thus, the 
predictions of the instanton model, under the considered condition a g = p c Vg < 1> are m qualitative agreement with 
information extracted from the lattice data. 



IV. CONCLUSIONS 



The instanton model provides a way of constructing the nonlocal vacuum condensates. We have obtained the 
expressions for the nonlocal gluon <: TrF^(x)E(x,y)F pa '(y)E(y,x) :> correlator beyond the single instanton (SI) 
approximation [^5|. They have consistent properties at short as well as at large distances. The model constructed 
predict the behaviour of nonperturbative part of gluon correlation functions in the short and intermediate region 
assuming that it is dominated by instanton vacuum component. 

To this goal, we have suggested that the instanton A^ 1 {x) is developed in the physical vacuum field b^x) inter- 
polating large-scale vacuum fluctuations. We have found that at small distances the instanton field dominates, and 
at large distances it decreases exponentially. We did not assume any particular properties of the long-wave vacuum 
field b^{x) but managed to reduce the effect to certain phenomenological quantities, namely, the correlation function 
B(x 2 ) determined by its strength (F£} b and the correlation length R. Within this model, by averaging over random 
color vector orientations of the background field with respect to the fixed instanton field orientation, we have found 
equation ( p0| ) governing the deformation of the instanton under the influence of the weak background vacuum field. 
Following Affleck idea we have shown that, to stabilize the instanton, we need to put constraints on the system. Next, 
we have found the constraint independent asymptotics of the instanton solution at large distances, given by Eqs. ( |29| ) 

and (|30|) , where it is exponentially suppressed A^ Ia {x) ~ 2rf^„x v — c - „ , exp — - (jj g \x\f'^ 2 unlike the powerful 

x 1 | x | L 3 J 

decreasing SI. It is important to note that the form of this asymptotics is also independent on the model for the 

( N c 2\ V 

background field and the driven parameter r/ g ~ I —R (F b J, I only weakly depends on it. Assuming that 

\y (iv c — 1) / 

the external field is weak, the CI profile function is close to SI profile at distances smaller than p c and it decreases ex- 
ponentially at distances larger than r]~ 1 (see (jl])). In particular, this result means that the leading interaction term of 
a widely separated instanton-anti-instanton pair in physical vacuum decays exponentially with separation and differs 
from dipole interaction term found previously in unconstrained model. The knowledge of the constraint-independent 
parts of CI allowed us to construct the solution in the ansatz form (|3^) with the profile functions ( p3| ) - ( f45| ) . As it is 
seen from Figs. (||) - (Q) the profile of the CI and its action are practically independent of the choice of the ansatz if 
the interference parameter is in the region, p c r\ g < 1, where our considerations are justified. 

Then, for an arbitrary classical gauge field of the form A^ Ia (x) = 2rj^ fl x u (p g (a; 2 ), we have found the expressions 
( p5| ) and ([56 ) for the combinations of form factors D (x 2 ) and D\ (x 2 ) , which parameterize the gauge-invariant gluon 
field strength correlator. These expressions generalize the previously known expressions for the SI model p5[. The 



14 



correlators have been calculated numerically. As it turns out, at a reasonable set of parameters, guaranteeing the 
smallness of the large-scale vacuum field fluctuations, the D (or) structure is close to the SI induced function with 
the exponential asymptotics being developed at large distances^]. At the same time, the D\ (x 2 ) structure is about 
two orders smaller than the D (x 2 ) function at any reasonable choice of the parameter p c r\g- As it is explained 
in the Appendix A, the reason is that in the dilute vacuum the CIs are "almost" self-dual. The relative strength 
of the form factors D (x 2 ) and D\ (x 2 ) is very sensitive to the accepted physical picture. The lattice data are 
in qualitative agreement with predictions of the constrained instanton model. It means, in particular, that in the 
interpretation of the lattice data more justified parameterization for the correlation functions can be used. It allows 
one to extract from data the values of physical interest and separate the perturbative tail from the nonperturbative 
contribution. Moreover, due to fast decay of the CI induced part of correlators, the exponential decay observed in 
lattice calculations can be attributed to the background component of the vacuum field, or be described by some other 
field theoretical approaches |]||. From the other side, the SI model is inconsistent with large distance behavior. The 
nonperturbative part of the functions A(x 2 ) and B{x 2 ) are the sum of short-range instanton induced contributions 
(|55|) and ((5^), multiplied by the weight factor n c 32ir 2 / (0 \F 2 \ OV ,, and the long-range contribution (|l9|), modeled 
by exponentially decreasing function (jH) with the weight factor (F 2 ") b / (0 |F 2 | OV . ,. 

The constrained instanton model introduces two characteristic scales (correlation lengths). One is related to short 
distance behavior of the correlation functions and another with long range distance behavior. The first one, A~ , 
is predictable and expressed in terms of physical quantities. In SI approximation, given by Eqs. (0) and (||), it is 
proportional to the instanton size, and gains small negative corrections due to the background, Eq. (j63|). As to large 
scale it is out of the scope of our model, and we can only physically relate it to the confinement size or extract it from 
long distance asymptotics of lattice calculations. The microscopic description of the long distance background field 
needs other considerations not examined in the present work. 

The calculations have been performed in a gauge - invariant manner by using the expressions for the instanton field 
in the axial gauge. The behavior of the correlation functions demonstrates that in the single constrained instanton 
approximation the model of nonlocal condensates can well reproduce the behavior of the functions at short and 
intermediate distances, while the large-scale asymptotics is dominated by the background field. It would be quite 
desirable to make a fit to the lattice data using as an input the instanton induced correlators. The important question 
concerning the interacting ensemble of the constrained instantons has also to be postponed for another specific work. 
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APPENDIX A: APPENDIX 

The question may arise why the short-range vacuum contributions to the D\ (x 2 ) structure are negligible. As it 
was found in || in the SI case D\ (x 2 ) = due to self-duality of instanton solution. The CIs are not self-dual and 
contribute to D\ (x 2 ), but to what extent the self-duality is violated? We are going to show that for the reasonable 
set of parameters and ansatze assumed the CIs are "almost" self-dual and this is the reason for smallness of D\ (x 2 ). 
On the contrary if lattice simulations detected very big contribution to D\ (x 2 ), it would mean that self-duality is 
lost and there is no chances to save instantons as individual objects in the QCD vacuum. From the point of view of 
our model, any essential contributions to D\ (x 2 ) can arise only from the large-scale vacuum fluctuations. 

The dual field strength F^J' a (x) = ^^v P aF^' a (x), where F^J> a (x) is defined in (|34|), can be expressed in the 
form 



5 The similar behavior is expected for the quark correlator in the physical vacuum Era 
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FpJ' a (s) = 4 [^,£3i (x) + (x^ p - a:^) x p cj 2 (x)] 



(Al) 



with forms u> (x) 



wi (re) = v s (z ) + i , ~z v - v - ~, v ~, - r sr y • fe2 



w 2 (a;) = ^2 (x) = (x 2 ) + 



Let us consider the difference of the field strengths given in the regular at zero gauge 



(A2) 



F% a {x)-F% a {x) 



% v + 2 



[XfiVvp x v 7 lp,p) Xp 



req i \ 2 
L0 2 (X) X 



The self-duality condition Ff?J a (x) — F^J a (x) = is satisfied for the SI case, where 



2, ,reg,I 



X UJ- 



Comparing Eqs. (34) and (Al) with Eq.([A3j) we can consider the condition 



(A3) 



(a:) = (see Eq. (|T|)) 



reg,CI 



(x) 



x 2 « 



reg,CI 



(x) 



CO 



reg,I 



(x) 



as a criterion indicating that the field is "almost" self-dual. It can be checked numerically, that it is really the case 
at reasonable choice of background field strength parameter a g < 1 and all forms of ansatze for CI. At the same time 
at larger values of parameter a g ~ 4 4- 6 the inequality is not fulfilled. Thus we show that assuming diluteness of the 
vacuum the CIs are "almost" self-dual solutions and, as a consequence, contribute very small to the D± (a; 2 ) structure 
of the vacuum gluon field strength correlators. 



APPENDIX B: APPENDIX 



The function <E>(z 2 ) for the different forms of the form factor B (a; 2 ): 



®g{z 2 ) = ^ [2V^V 3 erf (y) - 3y 2 + 1 - (l - 2y 2 ) exp (~y 2 )] , Gaussian, 

y 



and 



®m(z 2 ) = —j [4y 3 arctan (y) + y 2 - (l + 3y 2 ) In (l + y 2 )] , monopole, 

y 



<f>E(z 2 ) = [2y 3 - 3y 2 + 6 - 6 (1 + y) exp (-y)] , exponential, 

where y = z / R, R being the correlation length of the large-scale vacuum field. 

For completeness we present here the small interference contribution to the functions A(x 2 ) and B(x 2 ) 



(Bl) 
(B2) 

(B3) 



A^(x 2 ) = N D (0\F 2 \0) 6{N ,,_ 1)Jo 
■ (z + ■ Z-) {(z + ■ Z-) (3 — 4sin 2 (o!2)) — irr sin(2a 2 )} , 



drr 2 / dt${z + , z_)ip g (z 2 ) if g [z 2 _ ) 



(B4) 



where 



B mterf {x 2 ) =N D {0\F 2 \0) 



oo />oo 
2 



6(JVf-lW ^ Jo * $ (*+'*-)M4)M*-) 



•r 2 {2 (z+ ■ z-) (2 - 3sin 2 (a 2 )) - xrem(2a z )} , 



<f>(z + ,z_)=4 da d(3a(iB (az + - /?z_) 

z± = (r, i± z = (r, i) , and tp g (z 2 ) is defined in ((3^). In deriving the above expressions we used the Schwinger 
- Fock gauge for background field, Eqs. (0) - (|l9|), and neglected the derivative D'^x 2 ) in comparison with the 
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function B (a; 2 ) itself. In fact, it is consistent with D\ (x 2 ) = that follows from the lattice data and instanton 
calculations pq| . We find that the interference contributions to the correlators are very small in absolute value, have 
shorter correlation length comparing with the CI contributions ( |55|) , ( |56| ) and do not lead to valuable appearance of 
the Di (x 2 ) structure. 
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FIGURE CAPTIONS: 



Fig. 1: The constrained instanton profile functions x 2 tp g {\x\ j p) (|32|), corresponding to the ansatz (|45|), at different 
values of the parameter {prig) 2 : {pr) g ) 2 — solid line (instanton case), (pt] g ) 2 — 0.5- short dashed line, {prj g ) 2 = 3 - 
long dashed line. 

Fig. 2: The instanton (solid line) and constrained instanton profile functions x 2 y g (|a;| / p) ( |32"| ) corresponding to 
different ansatze: (^) dotted line, ([lij) - short dashed line, ([l5|) - long dashed line, at the large value of the parameter 

(pv 3 ) 2 ■■'>■ 

Fig. 3: The density action of the instanton (solid line) and constrained instanton G 2 (x 2 ) = p 4 /192F^ (x) F*£ (x), 
Eq. (p7|), as a function of |z| / p corresponding to different ansatze: ( (43f ) dotted line, (Q) - short dashed line, p^ ) - 
long dashed line, at the value of the parameter (prj g ) 2 = 1. 

Fig. 4: The classical action of the instanton (solid line) and the constrained instantonSc/, Eq. (|36|), as a function 
of prj g corresponding to different ansatze: d43| ) dotted line, (Q) - short dashed line, ([|5]) - long dashed line. The 
action is given in units of 8TT 2 /g 2 . 

Fig. 5: The form factors D (top lines) and D\ (bottom lines) (all normalized by D(0)) versus physical distance x, 
for the instanton size p = 0.3 fm and parameters (pi] g ) 2 — (solid lines) and (pr] g ) 2 = 1 (dashed lines). 

Fig. 6: The form factor D(p) as a function of p\p\ corresponding to the ansatz ((45|), at different values of the 
parameter (pf] g ) 2 : {prig) 2 — solid line (instanton case), [prj g ) 2 = 0.5 - short dashed line, {prj g ) 2 = 3 - long dashed 
line. 

Fig. 7: The interference term contribution to the gluon condensate normalized by the instanton contribution 
F^i/fA^j-Fpo-lyl 7 ] = n c 327r 2 as function of the large-scale vacuum fluctuation correlation length 1/(3 = R/p c and its 
strength parameter a g — i] g p c . 
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Fig. 1. 



FIG. 1. The constrained instanton profile functions x 2 ip g (\x\ / p) (p2|), corresponding to the ansatz (piq), at different values 
of the parameter {prig) 2 : (pr/ g ) 2 = solid line (instanton case), (pr/g) = 0.5 - short dashed line, (prig) 2 = 3 - long dashed line. 
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FIG. 2. The instanton (solid line) and constrained instanton profile functions x 2 (p g (\x\ / p) (p2[) corresponding to different 



ansatze: ( fl3| ) dotted line, ([5]) - short dashed line, (^|) - long dashed line, at the large value of the parameter {prig) 2 = 3 
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FIG. 3. The density action of the instanton (solid line) and constrained instanton G 2 (x 2 ) = p 4 /192_ F/f° ( x) Ff,?, (x), Eq. 
(^), as a function of \x\ / p corresponding to different ansatze: ( |43| ) dotted line, ( ji^ ) - short dashed line, ( jig ) - long dashed 
line, at the value of the parameter (prj g ) 2 = 1. 
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FIG. 4. The classical action of the instanton (solid line) and constrained instantonSc/ , Eq. (pq), as a function of pry 9 
corresponding to different ansatze: (^) dotted line, (^) - short dashed line, (^) - long dashed line. The action is given in 
units of 8n 2 /g 2 . 
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FIG. 5. The form factors D (top lines) and D\ (bottom lines) (all normalized by -D(O)) versus physical distance x, for the 
instanton size p = 0.3 fm and parameters {prj g ) 2 = (solid lines) and (pr/ g ) 2 = 1 (dashed lines). 
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FIG. 7. The interference term contribution to the gluon condensate normalized by the instanton contribution 
i ? (JI ,[j4 / ]Fp (T [yl / ] = n c 32-K 2 as function of the large-scale vacuum fluctuation correlation length 1//3 = R/ p c and its strength 
parameter a g — r\ g p c . 
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